Dimer mean-field model for the Ising spin glass 
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A dimer mean-field model for the Ising spin-glass is presented. Despite its simplicity it captures 
some of the essential features of the spin-glass physics. The distribution of the single-spin mag- 
netization is determined from a self-consistent integral equation. By solving the self-consistency 
condition numerically, we find that there are two temperature scales characterizing the glass tran- 
sition. At the first, higher temperature, the glass order parameter becomes non-vanishing, and at 
the second, freezing temperature, it saturates to its maximal value. The effect of magnetic field and 
the existence of the Almeida-Thouless line are discussed. Finally, it is shown that the information 
compressibility, defined as the derivative of entropy with respect to energy, diverges at the freezing 
temperature. This indicates a zero internal temperature and true glassy dynamics with diverging 
relaxation times. 

PACS numbers: 75.10.Nr,75.50.Lk 
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The interest in spin-glasses (SGs), by now several 
decades old 1-4 , stems from the fact that in spite of their 
seemingly simple description, they display a rich and 
complex behavior 5 . Specifically, quantum SGs have been 
investigated using elaborate mathematical formalisms 
such as the replica approach 3 and quantum annealing 
methods 6 . Studying the exact properties from exact 
diagonalization of the Hamiltonian requires tremendous 
computational effort, since the dimension of the Hilbcrt 
space increases exponentially with system size 7,8 . There- 
fore, it has long been recognized that simplified mod- 
els, which exhibit the glass properties but are easily ad- 
dressed and solved are highly desirable 9-11 . 

Here, we present such a model based on the short-range 
Ising spin-glass model. We study its single-spin proper- 
ties, using a dimer-mean-field approximation for the dis- 
tribution of the single-spin magnetization. The model is 
simple and physically transparent, and yet captures some 
of the essentials of the spin-glass transition. Our main 
results are as follows, (i) We find that in D > 2 dimen- 
sions there are two important temperature scales. The 
first, higher temperature, corresponds to the onset of a 
finite order parameter. The second temperature, is one 
at which the order parameter saturates to its maximal 
value. We can thus distinct between the "glass transi- 
tion" , at which the order parameter first appears, and 
the "freezing transition" , at which it saturates. Inter- 



estingly, we find that in ID there is only a glass tran- 
sition but no freezing transition, (ii) For D > 2, there 
exists a region of magnetic fields in which the freezing 
transition is supported, thus confirming the existence of 
an Almeida-Thouless line 3 ' 5 ' 8 ' 12 for this model, (iii) We 
show that the internal temperature - as experienced by 
the spins - vanishes at the freezing transition. We argue 
that consequently the dynamics exhibits a critical slow- 
ing down. 

Model and method of calculation - The starting point 
of our calculation is the Anderson-Edwards Hamiltonian 
for the short-range Ising spin glass 1 , 
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where (ij) indicates sum over nearest neighbors on a 
square D-dimensional lattice, are the values ± J ran- 
domly given with equal probability, and Sj is the z- 
direction spin operator for the particle at position i. 

The next step is to apply a dimer-mean-field approxi- 
mation. We divide the system into pairs of spins (spin- 
dimers, see Fig. 1), and we write down explicitly the 
Hamiltonian of the spin-dimer. In the basis of the dimer 
states | IT), I Tl))| IT},! II)) the Hamiltonian is simply 
given by 
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where hi and h 2 are the magnetic fields experienced by J12 is the coupling of the latter. The mean-field approx- 
the first and second spins of the dimer, respectively, and imation is employed by setting hi = h ex ^ — Ji(z — l)sj, 
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where ft ex t is the external magnetic field (in the z- 
direction), J 4 is a (random) exchange coupling, z is the 
number of nearest neighbors, and Sj is the (average) spin 
neighboring the spin s, (see Fig. 1). In fact, for a more 

exact treatment one should write hi = ^ext+X/fci ^ij^ji 
that would amount to accounting for each of the neigh- 
boring spins separately. However, we found by comparing 
numerical calculations on both versions (we used a two- 
dimensional square lattice for the second) that this does 
not alter the results. 
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FIG. 1: (Color online) A schematic of the dimer mean- field 
model. Out of the whole spin lattice, two neighboring spins 
are chosen, and their Hamiltonian is written, where the ex- 
ternal field experienced by the spins results from the average 
field generated by their neighboring spins. 

Since the Hamiltonian is diagonal, it is now a matter 
of straightforward algebra to evaluate the average spin 
of either member of the spin-dimer, which is given by 
(setting kg = ft = 1) 
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where Z = Tr(e~ w/T ) = X« =1 e £l/T with the diago- 
nal elements of the Hamiltonian, and T is the external 
temperature. 

Since all the Js in the above expression are random 
variables, s and s are also random variables. Therefore, 
it is the distribution of s which should be determined sclf- 
consistently. Since the function s(J i2 , J\, J 2 , s\, s\, h ext ) 
is known, the equation for the distribution function of s, 
g(s) reads 

g(s) = J dJi 2 dJ 1 dJ 2 di 1 dS 2 /(J 12 )/(J 1 )/(J 2 ) x 

'xg(si)g{s 2 )S{s - s( J12, Ji, J 2 , h,h,h eXt )) (4) 

where S is the Dirac delta-function, and f(Ji) = 
\ (S(Ji — J) + S(Ji + J)) is the distribution function for 
the couplings. 

Eq. 4 is a central result of this paper and it represents 
the mean-field self-consistent equation for this model. In 
the following, we solve it numerically for different tem- 
peratures, external fields, etc. From the solution of g(s), 
we can estimate all single-spin quantities, such as the av- 
erage magnetization (s), the so-called Edwards- Anderson 
glass order parameter, which in this mean-field approx- 
imation is simply (s 2 ) and other properties such as the 



entropy, energy, etc. Here, we note that naively one could 
suggest a similar method for a single spin instead of a 
dimer, yet it turns out that this cannot capture the basic 
spin- glass properties. The reason is that the glassy be- 
havior is a direct consequence of the competition between 
disorder and interactions, and thus a dimer structure is 
the minimal one that includes both. 

The glass and freezing transitions - We start by cal- 
culating the distribution function g(s) for different tem- 
peratures T (scaled to units of J), with no external field, 
for a two-dimensional (2D) system. In Fig. 2 the distri- 
bution function is plotted for temperatures T/J = 1, 2, 3 
and 3.9. For T/J > 4 (which is the number of nearest 
neighbors, z, in 2D) the distribution becomes strongly 
peaked around s = (not shown). 




FIG. 2: (Color online) Distribution function g(s) for a two- 
dimensional spin glass for temperatures T/J — 1, 2, 3 and 3.9. 

As seen from Fig. 2, the distributions are symmetric 
around s — for all temperatures (up to numerical er- 
rors), and hence (s)(t) — 0- More information about 
the single spin may be obtained by evaluating the second 
moment, (s 2 ), the spin-glass order parameter. 

In Fig. 3 we plot (s 2 ) as a function of T/J for the 2D 
spin glass. One can observe two relevant temperature 
scales. The first, Tq = zJ, corresponds to the onset 
of a non-zero (s 2 ), and we thus recognize this point as 
the "glass transition". However, at T{ = J the order 
parameter reaches (s 2 ) — 1, which implies that thermal 
fluctuations are completely suppressed. This point thus 
corresponds to a "freezing transition" . 

In the inset of Fig. 3 we plot the same order param- 
eter for 3D and ID glasses (note that in this mean-field 
model differences originate only from the number of near- 
est neighbors z). The temperature is now scaled by T/zJ, 
and the arrows indicate the point at which T/J =1. In 
3D the behavior is similar to that of 2D, but in ID we find 
that there is no freezing transition, and (s 2 ) remains be- 
low (s 2 ) = 1 down to the lowest temperatures reachable 
numerically, indicating that the system never "freezes" . 
Thus, for this model the critical dimension for the freez- 
ing transition is D = 2. 

Effect of magnetic field - A standing question in the 
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FIG. 3: (Color online) (s 2 ) as a function of T/J for the 2D 
spin glass. At a temperature Tg = zj, (s 2 ) becomes nonzero, 
indicating a glass transition. At Tt = J, (s 2 ) ~ 1 indicating 
the suppression of thermal fluctuations and a freezing tran- 
sition. Inset: the same for 3D and ID. Arrows indicate the 
position of Tf — J. As seen, in 3D the behavior is similar to 
that of 2D, but in ID there is no freezing transition. 



theory of spin-glasses is the effect of an external mag- 
netic field 8 . In a seminal paper, Almeida and Thouless 
predicted that the (mean-field) spin-glass phase is ro- 
bust against small external fields, and hence the system 
still has a glass phase and exhibits a glass transition 12 . 
Large fields, however, polarize the spins regardless of the 
random interactions, and the glassy behavior vanishes. 
Thus, there is a crossover from glassy to polarized ground 
state, the so-called Almeida-Thouless line. It is of in- 
terest to see whether our model displays these features, 
which are supposed to be common to all mean-field SGs. 

In Fig. 4 we show the single-spin magnetization (s) in 
two dimensions, evaluated from the self-consistent solu- 
tion of Eq. 4 for finite temperature and magnetic field h 
(measured in units of J). Fig. 4(a) shows (s) as a func- 
tion of temperature for various values of h. For small 
h, as the temperature approaches the freezing tempera- 
ture, the glassy behavior begins to dominate and (s) is 
suppressed. For larger values of h, we find finite magne- 
tization even at the freezing transition. In Fig. 4(b) we 
plot (s) at the freezing transition. As seen, at around 
h ~ J the spins develop finite magnetization, indicating 
the existence of an Almeida-Thouless line for this model. 

Effective temperature and diverging relaxation times - 
In this final section we wish to study the relation between 
the entropy S and energy E for the present model. As 
we will show, their relation leads to a vanishing internal 
temperature at the freezing transition and we will argue 
that this implies diverging relaxation times, which are 
one of the hallmarks of glassy behavior. 

In a recent paper 13 , we have introduced the dynam- 
ical quantity Ki(t) = ^f(^f)" 1 |t'=t (we set k B = l), 
which we call the information compressibility. We have 
shown that for a two-level system open to environments, 
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FIG. 4: (Color online) (a) The single-spin magnetization (s) 
(in two dimensions) as a function of temperature for different 
values of external magnetic field, (b) The single-spin magne- 
tization (s) at T/J — 1 as a function of external fields. At 
small fields the magnetization vanishes, and becomes finite at 
about h ~ J, indicating the existence of an Almeida-Thouless 
line. 

this quantity defines a dynamic temperature in non- 
equilibrium processes, by the simple relation T e g = kJ 
(at equilibrium, T c g coincides with the thermodynamic 
temperature), as the temperature of a probe adjusted 
so that the system's properties are minimally perturbed. 
When possible, one can use the chain rule and define k/ 
also for steady-states as Kj = dS/dE. This latter quan- 
tity is well studied in simple glass models and is known 
to vanish at the transition 9,11 . 

In our model, both the energy and the entropy are eas- 
ily calculated from the Hamiltonian (Eq. 2). Of course, 
they are both random variables and depend on the spe- 
cific realizations of the neighboring spins and the cou- 
pling constants. We thus calculate the average values of 
S and E as a function of temperature (at zero magnetic 
field) , with the use of the distributions obtained earlier 
(Fig. 2). From these we can construct the function S(E). 

In Fig. 5 we plot T cff = iff 1 = (dS/dE)' 1 as a func- 
tion of temperature T. As can be seen, the effective 
temperature vanishes at the freezing transition (T = J). 
This provides an intuitive explanation for the suppres- 
sion of thermal fluctuations: a small change in energy 
induces a large change in entropy, the number of mi- 
crostates available to the system increases, and therefore, 
the system can explore an ever increasing number of mi- 
crostates thus reducing the relative fluctuations. This is 
equivalent to stating that close to the freezing transition 
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fluctuations are dominated by the effective temperature 
T e g (which vanishes at the transition) rather than the 
external temperature. 
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FIG. 5: The effective temperature, T e « — (dS/dE) -1 , as a 
function of T (in two dimensions). The effective tempera- 
ture vanishes at the freezing transition, giving an intuitive 
explanation for the suppression of thermal fluctuations and 
diverging relaxation time at that point (see text). 

To understand how a vanishing T c g leads to a diverging 
relaxation time (defined as the time it takes the system 
to restore its equilibrium state after some excitation) we 
turn to a model introduced by Majumdar et ai. 14 These 
authors have introduced a general model for glassy dy- 
namics, in which the system is described by an order pa- 
rameter p (in our case, that corresponds to (s 2 )). They 
then considered a situation in which transition rates be- 
tween states of different values of p (say p and p') are 
determined not by the usual Boltzmann rates, but rather 
by the following rule: 

W P-f-\ e -AE/T j AE>Q (5) 

where AE = E(p') - E[p), and AS = S{p) - S(p'). The 
authors of Ref. 14 have showed that considering such 
dynamics leads directly to diverging relaxation times 
at the transition, since the number of phase-space tra- 
jectories that correspond to a decrease in energy de- 
creases as one approaches the glass (or freezing) transi- 



tion. To relate to our finding, we note that for processes 
which involve a small energy difference, one can write 
AS" ~ j^AE = AE/T eS . This means that, due to the 
interplay between the change in temperature and change 
in the number of available microstates, the temperature 
which defines the relaxation processes - that relax higher 
energy states into lower ones - essential for equilibrium 
restoration is T e g rather then the external temperature. 
As T e g — > 0, the rates to go downhill in energy vanish, 
thus resulting in diverging relaxation times. 

Summary. - In summary, we have presented a sim- 
ple dimer mean-field model to calculate the single-spin 
properties in the Ising spin-glass. This model, which is 
physically transparent and computationally cheap, gives 
insight into the properties of the glass transition. Specif- 
ically, we have shown that at the freezing temperature 
Tf= J thermal fluctuations are suppressed. We have also 
demonstrated the effect of magnetic field, and showed 
that the effective temperature of the spins vanishes at the 
freezing transition. This, in turn, provides an intuitive 
explanation of the suppression of thermal fluctuations 
and divergence of timescales at the freezing transition of 
the Ising spin-glass. 

In the future we plan to extend the model in various 
directions. One possible direction is to include more than 
nearest neighbors, and by using real-space renormaliza- 
tion methods we plan to obtain information not only 
about the single spin but also about correlations between 
distant spins, which are very important in characteriz- 
ing the glassy state. Other directions include extending 
the system to three-dimensional spins (e.g., Heisenberg 
spin-glasses) , including magnetic fields in a general direc- 
tion, and extending the calculation to larger spins (rather 
than spin i). These studies can be relevant to various 
systems such as memory-shape alloys 15,16 and strained 
materials 17 . 

We thank F. Krzakala for valuable discussions. This 
research was funded by DOE under grant DE-FG02- 
05ER46204. 
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